L\'evy flights and self-similar exploratory behaviour of termite workers: beyond model fitting

Animal movements have been related to optimal foraging strategies where self-similar trajectories are central. Most of the experimental studies done so far have focused mainly on fitting statistical models to data in order to test for movement patterns described by power-laws. Here we show by analyzing over half a million movement displacements that isolated termite workers actually exhibit a range of very interesting dynamical properties --including L\'evy flights-- in their exploratory behaviour. Going beyond the current trend of statistical model fitting alone, our study analyses anomalous diffusion and structure functions to estimate values of the scaling exponents describing displacement statistics. We evince the fractal nature of the movement patterns and show how the scaling exponents describing termite space exploration intriguingly comply with mathematical relations found in the physics of transport phenomena. By doing this, we rescue a rich variety of physical and biological phenomenology that can be potentially important and meaningful for the study of complex animal behavior and, in particular, for the study of how patterns of exploratory behaviour of individual social insects may impact not only their feeding demands but also nestmate encounter patterns and, hence, their dynamics at the social scale.Comment: 13 pages, 11 figures. Unrevised version. Final version to appear in Plos ON

A quantitative model of human jejunal smooth muscle cell electrophysiology

10.1371/journal.pone.0042385PLoS ONE78

Normalized Lift: An Energy Interpretation of the Lift Coefficient Simplifies Comparisons of the Lifting Ability of Rotating and Flapping Surfaces

For a century, researchers have used the standard lift coefficient CL to evaluate the lift, L, generated by fixed wings over an area S against dynamic pressure, ½ρv2, where v is the effective velocity of the wing. Because the lift coefficient was developed initially for fixed wings in steady flow, its application to other lifting systems requires either simplifying assumptions or complex adjustments as is the case for flapping wings and rotating cylinders

Wake Development behind Paired Wings with Tip and Root Trailing Vortices: Consequences for Animal Flight Force Estimates

Recent experiments on flapping flight in animals have shown that a variety of unrelated species shed a wake behind left and right wings consisting of both tip and root vortices. Here we present an investigation using Particle Image Velocimetry (PIV) of the behaviour and interaction of trailing vortices shed by paired, fixed wings that simplify and mimic the wake of a flying animal with a non-lifting body. We measured flow velocities at five positions downstream of two adjacent NACA 0012 aerofoils and systematically varied aspect ratio, the gap between the wings (corresponding to the width of a non-lifting body), angle of attack, and the Reynolds number. The range of aspect ratios and Reynolds number where chosen to be relevant to natural fliers and swimmers, and insect flight in particular. We show that the wake behind the paired wings deformed as a consequence of the induced flow distribution such that the wingtip vortices convected downwards while the root vortices twist around each other. Vortex interaction and wake deformation became more pronounced further downstream of the wing, so the positioning of PIV measurement planes in experiments on flying animals has an important effect on subsequent force estimates due to rotating induced flow vectors. Wake deformation was most severe behind wings with lower aspect ratios and when the distance between the wings was small, suggesting that animals that match this description constitute high-risk groups in terms of measurement error. Our results, therefore, have significant implications for experimental design where wake measurements are used to estimate forces generated in animal flight. In particular, the downstream distance of the measurement plane should be minimised, notwithstanding the animal welfare constraints when measuring the wake behind flying animals

Constructing a Stochastic Model of Bumblebee Flights from Experimental Data

PMCID: PMC3592844This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Comparing stochastic differential equations and agent-based modelling and simulation for early-stage cancer

There is great potential to be explored regarding the use of agent-based modelling and simulation as an alternative paradigm to investigate early-stage cancer interactions with the immune system. It does not suffer from some limitations of ordinary differential equation models, such as the lack of stochasticity, representation of individual behaviours rather than aggregates and individual memory. In this paper we investigate the potential contribution of agent-based modelling and simulation when contrasted with stochastic versions of ODE models using early-stage cancer examples. We seek answers to the following questions: (1) Does this new stochastic formulation produce similar results to the agent-based version? (2) Can these methods be used interchangeably? (3) Do agent-based models outcomes reveal any benefit when compared to the Gillespie results? To answer these research questions we investigate three well-established mathematical models describing interactions between tumour cells and immune elements. These case studies were re-conceptualised under an agent-based perspective and also converted to the Gillespie algorithm formulation. Our interest in this work, therefore, is to establish a methodological discussion regarding the usability of different simulation approaches, rather than provide further biological insights into the investigated case studies. Our results show that it is possible to obtain equivalent models that implement the same mechanisms; however, the incapacity of the Gillespie algorithm to retain individual memory of past events affects the similarity of some results. Furthermore, the emergent behaviour of ABMS produces extra patters of behaviour in the system, which was not obtained by the Gillespie algorithm

Diving-flight aerodynamics of a peregrine falcon (Falco peregrinus)

This study investigates the aerodynamics of the falcon Falco peregrinus while diving. During a dive peregrines can reach velocities of more than 320 km h⁻¹. Unfortunately, in freely roaming falcons, these high velocities prohibit a precise determination of flight parameters such as velocity and acceleration as well as body shape and wing contour. Therefore, individual F. peregrinus were trained to dive in front of a vertical dam with a height of 60 m. The presence of a well-defined background allowed us to reconstruct the flight path and the body shape of the falcon during certain flight phases. Flight trajectories were obtained with a stereo high-speed camera system. In addition, body images of the falcon were taken from two perspectives with a high-resolution digital camera. The dam allowed us to match the high-resolution images obtained from the digital camera with the corresponding images taken with the high-speed cameras. Using these data we built a life-size model of F. peregrinus and used it to measure the drag and lift forces in a wind-tunnel. We compared these forces acting on the model with the data obtained from the 3-D flight path trajectory of the diving F. peregrinus. Visualizations of the flow in the wind-tunnel uncovered details of the flow structure around the falcon's body, which suggests local regions with separation of flow. High-resolution pictures of the diving peregrine indicate that feathers pop-up in the equivalent regions, where flow separation in the model falcon occurred

Geometric Mixing, Peristalsis, and the Geometric Phase of the Stomach

Mixing fluid in a container at low Reynolds number - in an inertialess environment - is not a trivial task. Reciprocating motions merely lead to cycles of mixing and unmixing, so continuous rotation, as used in many technological applications, would appear to be necessary. However, there is another solution: movement of the walls in a cyclical fashion to introduce a geometric phase. We show using journal-bearing flow as a model that such geometric mixing is a general tool for using deformable boundaries that return to the same position to mix fluid at low Reynolds number. We then simulate a biological example: we show that mixing in the stomach functions because of the "belly phase": peristaltic movement of the walls in a cyclical fashion introduces a geometric phase that avoids unmixing.Comment: Revised, published versio

On the Mysterious Propulsion of Synechococcus

We propose a model for the self-propulsion of the marine bacterium Synechococcus utilizing a continuous looped helical track analogous to that found in Myxobacteria [1]. In our model cargo-carrying protein motors, driven by proton-motive force, move along a continuous looped helical track. The movement of the cargo creates surface distortions in the form of small amplitude traveling ridges along the S-layer above the helical track. The resulting fluid motion adjacent to the helical ribbon provides the propulsive thrust. A variation on the helical rotor model of [1] allows the motors to be anchored to the peptidoglycan layer, where they drive rotation of the track creating traveling helical waves along the S-layer. We derive expressions relating the swimming speed to the amplitude, wavelength, and velocity of the surface waves induced by the helical rotor, and show that they fall in reasonable ranges to explain the velocity and rotation rate of swimming Synechococcus

Effects of X-ray dose on rhizosphere studies using X-ray computed tomography

X-ray Computed Tomography (CT) is a non-destructive imaging technique originally designed for diagnostic medicine, which was adopted for rhizosphere and soil science applications in the early 1980s. X-ray CT enables researchers to simultaneously visualise and quantify the heterogeneous soil matrix of mineral grains, organic matter, air-filled pores and water-filled pores. Additionally, X-ray CT allows visualisation of plant roots in situ without the need for traditional invasive methods such as root washing. However, one routinely unreported aspect of X-ray CT is the potential effect of X-ray dose on the soil-borne microorganisms and plants in rhizosphere investigations. Here we aimed to i) highlight the need for more consistent reporting of X-ray CT parameters for dose to sample, ii) to provide an overview of previously reported impacts of X-rays on soil microorganisms and plant roots and iii) present new data investigating the response of plant roots and microbial communities to X-ray exposure. Fewer than 5% of the 126 publications included in the literature review contained sufficient information to calculate dose and only 2.4% of the publications explicitly state an estimate of dose received by each sample. We conducted a study involving rice roots growing in soil, observing no significant difference between the numbers of root tips, root volume and total root length in scanned versus unscanned samples. In parallel, a soil microbe experiment scanning samples over a total of 24 weeks observed no significant difference between the scanned and unscanned microbial biomass values. We conclude from the literature review and our own experiments that X-ray CT does not impact plant growth or soil microbial populations when employing a low level of dose (<30 Gy). However, the call for higher throughput X-ray CT means that doses that biological samples receive are likely to increase and thus should be closely monitored